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The reversion-of-series method is extended to the 
s - domain by utilizing non-linear Laplace transforms. 

The reversion of series in the s - domain is applied to a 
non-linear differential equation and approximate solutions 
are obtained. The approximate solution is modified for 
the case where the steady state is a constant value by 
calculating the exact steady-state value and applying it 
to the reversion approximation. The non-linear differen- 
tial equation considered is Duffing's equation with a 
damping term and sinusoidal and constant forcing functions. 
The theoretical solutions are compared to machine solutions. 
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I . INTRODUCTION 



Utilization of the Laplace transformation to obtain an 
accurate analytical solution of a non-linear differential 
equation has been prevented, in past, by the inability to 
obtain an appropriate expression for the Laplace transform 
of the non-linear terms. Baycura [Ref. 1] developed an 
expression for the Laplace transform of certain non-linear 
terms and Brady [Ref. 2] obtained a formulation for the 
general case. The non-linear transform is used in con- 
junction with the reversion-of-series method to obtain an 
approximate solution for a non-linear differential 
equation . 

The approximate solution is formulated and examined 
in detail by comparison to a machine solution. A solution 
obtained by modifying the approximate solution is also 
derived and examined in detail. 
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II. TRANSFORM EXPRESSION FOR INTEGRAL POWERS OF A FUNCTION 



The general formula for the Laplace transform of a 
function raised to an integral power [Ref. 3] is: 

I Jl tx n (t)l = s n- 1 X n (s) ( 

where n is a positive integer greater than zero. 



( 1 ) 
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III. REVERSION OF SERIES 



If a series is represented by 



2 ^ Z| 5 

x = a x y + a 2 y + a^y + a^y + a^y^ + ... (a^ 0) 



( 2 ) 



the coefficients of the series 



2 3 4 5 

y = A^x + AgX + A_x + A^x + A^x +... 



(3) 



are given by [Ref. 4] : 

A 1 = a^ A 2 = " ^3 A 3 = (2a 2 2 " a l a 3 5 



A rj ( 5a-2_ a 2 a 3 a l a 2j 5 a 2 ) 



A^ = — ^ (6a 1 2 a 2 a^ + 3a 1 2 a^ 2 + l4a 2 ^ - a^a^ - 21a^a 2 2 a2) 
a l 

Utilizing the relationship expressed in Equation (1) the 
Laplace transform of Equation (2) becomes: 



X(s) = a^Y(s) +a 2 sY 2 (s)+a^s 2 Y^(s)+a i |S^Y it (s) + ... 



(4) 



Let : 



2 3 

bf = a 1 b 2 = a 2 s b^ = a^s b^ = a^s 



Equation (4) can then be rewritten as: 



X(s) = b-^Y(s) + b 2 Y 2 (s) + b^Y^(s) + b^Y^(s) 



(5) 
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The Laplace transform of Equation (3) is: 

Y(s) = AjX ( s ) + A 2 sX 2 (s) + A 3 s 2 X 3 (s) + A 4 s 3 X 4 (s) ... (6) 

Let : 



2 2 

Bf = a-^ B 2 = ^2 S ~ ^3 S = ^4 S ' 

Equation (6) can then be rewritten as: 

Y ( s ) = B ^X ( s ) + B 2 X 2 (s) + B 3 X 3 (s) + B^X^ ( s ) ... (7) 



The coefficients of Equation (7) can be expressed in terms 
of the coefficients of Equation (5): 



B 1 ■ ^ B 2 - - B 3 = ^5 (2b 2 2 - b lV 



B 4 = (5b 1 b 2 b 3 - b x 2 b 4 - 5b 2 3 ) 

b l 



B 5 = — — q- (6b 1 2 b 2 b i| + 3b x 2 b 3 2 + l4b 2 4 - b x 3 b 5 - 21b 1 b 2 2 b 3 ) 



Thus, Y(s) can be expressed as a function of X(s) by the 
reversion of Equation (5) into Equation (7). 
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IV. GENERAL SOLUTION OF A NON-LINEAR DIFFERENTIAL EQUATION 



The technique illustrated in the previous section will 
now be used to solve Duffing 's equation with a damping 
term and initial conditions. The equation is expressed as: 



+ a it + w o y + hy3 = 



( 8 ) 



2 

where a, and h are constants and x(t) is a forcing func- 
tion. The Laplace transform of Equation (8) is: 



(s 2 + as + w 2 )Y(s) + hs 2 Y 3 (s) = X(s) + sy(0) + ||-| +ay(0) (9) 



Let : 



. 2 , , 2 
b^ = s + as + Wq 

^3 = 

Z (s ) = X(s) + y ( 0 ) ( s + a) + ^r| 



Equation (9) can then be written as: 

b^Cs) + b 3 Y 3 (s) = Z ( s ) (10) 

Consider Equation (10) as series with b n equal to zero for 
all n except n equal to one and n equal to three . 

By reversion of series: 



Y(s) 



Z ( s ) 

, 2 ^ ^ 2 , 

(s +as + coQ) 



hs 2 Z 3 (s) 3h 2 s i< Z 5 (s) 

p p IT + p p 7 

(s +as + ^Q) (s +os+Wq) 



(ID 
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The denominator of each term in Equation (11) can be 
factored. Thus: 



2 2 
s + as + (o q = 



s + 



a 

2 



2 2 2 ‘ 
(a-4co 2 ) 



i* _ - 



s + 



a 



( 2 u 2\ 2 ' 
(a -4 u)q) 



( 12 ) 



Let : 



A 





-4w 



2 

0 



) 



(13) 



B 





(14) 



Equation (11) can then be rewritten as: 

Yf'-'i - Z(s) hs 2 Z 3 (s) 3h 2 s i| Z 5 (s) , 

Y(s ) n w + 7 7 + . . . ( 

(s+A)(s+B) (s+A; (s+B; (s+A) (s+B) 

Equation (15) represents the general solution in the 
s - domain. Inspection of Equation (15) reveals that the 
form of the time solution depends on the nature of the 
forcing function and the relationship between a and 4(0 q . 
Further, if the forcing function is Laplace transformable 
and if h is equal to zero, the solution is exact. 
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V. APPLICATION OP THE GENERAL SOLUTION 



The general solution just derived will now be applied 
for zero initial conditions with constant and sinusoidal 
forcing functions . The time solutions are given in terms 
of numbered coefficients. The formula for each coefficient 
is listed in the Appendix. 

A. CONSTANT FORCING FUNCTION 

The inverse Laplace transform of Equation (15) yields 
three different solutions depending on the nature of the 
radical in Equations (13) and (14). 

1 . Radical Equal to Zero 

For the radical to be equal to zero, then: 




and 

A = B = | (17) 

The general solution in the s - domain is then: 

Y(s ) = 2 — hc3 » + 3h2c ^ + . . . (18) 

s(s+A) s(s+A) s(s+A) 

where c is the constant forcing function. 

The inverse Laplace transform of the first three terms of 
Equation (18) yields: 
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y ( t ) = K 1 + (K 2 +K 3 t)e At + + (K^ + Kgt+I^t 2 + . . . + K lQ t 13 )e At 

(19) 

where the , i = 1, 2, 3, . .., are defined In the Appendix. 

An examination of Equation (19) indicates that the initial 
conditions are satisfied. The effect of the various param- 
eters on the solution can be seen by examining the steady- 
state solution which can be expressed as: 
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It should be noted that only the first three terms of 
Equation (20) are utilized in Equation (19). Since Equation 
(20) is an alternating infinite series, a condition for 
convergence [Ref. 51 is that the magnitude of each successive 
term must decrease. This implies that the applicability 

of the reversion method is limited in the steady state to 

2 2 3 

the situation where he is small compared to (o)q) . 

2 

Using the arbitrary values c = 1.0, a = 4.0, oJq = 4.0 
and h = 1.0, the solution represented by Equation (19) 
is compared in Table I to a machine solution which was 
obtained using Runge-Kutta Adams - Moulton with error check. 
As the results indicate, a relatively accurate approximate 
solution was obtained. In Table II the solutions are com- 
pared for the coefficient of the non-linear term equal to 
ten. In this case the solution was a good approximation, 
but not as accurate as in the previous example. Examination 
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TABLE I 





a=4 . 0 , 


Wq= 4 . 0 , h=l . 0 , 


O 

rH 

II 

O 




Time 


Linear 


Reversion 


Machine 


Error 


0.0 


0.0000 


0.0000 


0.0000 


0.0000 


0.5 


0.0661 


0.0661 


0.0661 


0.0000 


1.0 


0.1485 


0.1485 


0.1484 


0.0001 


1.5 


0.2002 


0.2002 


0.1997 


0.0005 


2.0 


0.2271 


0.2269 


0.2259 


0.0010 


2.5 


0.2399 


0.2394 


0.2379 


0.0015 


3.0 


0.2457 


0.2447 


0.2430 


0.0017 


3.5 


0.2482 


0.2466 


0.2451 


0.0015 


4.0 


0.2492 


0.2471 


0.2458 


0.0013 


4.5 


0.2497 


0.2471 


0.2461 


0.0010 


5.0 


0.2499 


0.2467 


0.2462 


0.0005 


5 . 5 


0.2499 


0.2466 


0.2463 


0.0003 


6.0 


0.2500 


0.2465 


0.2463 


0.0002 


6.5 


0.2500 


0.2464 


0.2463 


0.0001 


7.0 


0.2500 


0.2463 


0.2463 


0.0000 


7.5 


0.2500 


0.2463 


0.2463 


0.0000 
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TABLE II 





a=4 . 0 j 


Wq= 4.0, h=10 . 0 , 


O 

1 — 1 
II 

O 




Time 


Linear 


Reversion 


Machine 


Error 


0.0 


0.0000 


0.0000 


0.0000 


0.0000 


0.5 


0.0661 


0.0661 


0 . 0660 


0.0001 


1.0 


0.1485 


0.1485 


0 . 1477 


0.0008 


1.5 


0.2002 


0.1997 


0.1956 


0.0031 


2.0 


0.2271 


0.2251 


0.2160 


0.0091 


2.5 


0.2399 


0.2347 


0.2221 


0.0126 


3.0 


0.2457 


0.2357 


0.2232 


0.0125 


3.5 


0.2482 


0.2327 


0.2229 


0.0098 


4.0 


0.2492 


0.2285 


0.2227 


0.0058 


4.5 


0.2497 


0.2246 


0.2225 


0.0021 


5.0 


0.2499 


0.2219 


0.2225 


-0.0006 


5.5 


0.2499 


0.2205 


0.2225 


-0.0020 


6.0 


0.2500 


0.2202 


0.2225 


-0.0023 


6.5 


0.2500 


0.2206 


0.2225 


-0.0019 


7.0 


0.2500 


0.2220 


0.2225 


-0.0005 


7.5 


0.2500 


0.2233 


0.2225 


0.0008 


8.0 


0.2500 


0 . 2246 


0.2225 


0.0021 


9.0 


0.2500 


0.2268 


0.2225 


0.0043 


15.0 


0.2500 


0.2292 


0.2225 


0.0067 
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of Equation (20) reveals that for h equal to ten, successive 
terms decrease in magnitude, but not as rapidly as in the 
case h was equal to one. If four terms of Equation (20) 
had been utilized in the solution, a steady-state value of 
0.2178 would have been obtained resulting in a better 
approximation. Table III represents a comparison of the 
solutions for the forcing function equal to ten. In this 
case the solution limits were exceeded and a large positive 
number was obtained for the steady state. If an additional 
term of Equation (20) had been included in the solution, a 
larger negative value would have been obtained for the 
steady state. Good results were realized however, for the 
initial portion of the transient solution, indicating the 
applicability of the reversion method in this region. The 
linear solution represented by the first two terms of 
Equation (19) was also included in the tables to illustrate 
a less accurate but simplified approximate solution. 

2 . Radical Greater than Zero 

2 2 

The relationship between a and 4 o>q is expressed as: 

a^>4cjQ . (21) 

Equation (15) is then expressed as: 



Y(s) = 



he 



ov 2 5 

3h c^ 



s(s+A )(s+B) 



s ( s+A) ^ ( s+B ) ^ s(s + A)7(s + B)7 



. ( 22 ) 



The inverse Laplace transform of the first two terms of 
Equation (22) yields: 
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TABLE III 





a=4 


.0, u)q = 4 . 0 j h-1.0. 


c=10 . 0 




Time 


Linear 


Reversion 


Machine 


Error 


0.0 


0.0000 


0.0000 


0.0000 


0.0000 


0.2 


0.1539 


0.1539 


0.1539 


0.0000 


0.4 


0.4780 


0.4780 


0.4770 


0.0010 


0.6 


0.8434 


0.8433 


0.8361 


0.0072 


0.8 


1.1877 


1.1867 


1.1628 


0.0239 


1.0 


1.4850 


1.4807 


1.4085 


0.0722 


1.5 


2.0020 


1.9560 


1.6365 


0.3195 


2 . 0 


2.2710 


2 . 0727 


1.5880 


0.4847 


2.5 


2.3990 


1.8907 


1.5006 


0.3901 


3.0 


2.4570 


1.5230 


1.5515 


-0.0285 


3.5 


2 . 48l8 


1.1488 


1.5567 


-0 . 4079 


4.0 


2.4925 


0 . 9814 


1.5575 


-0.5761 


4.5 


2.4969 


1.2081 


1.5569 


-0.3488 


5.0 


2.4988 


1.9344 


1.5567 


0.3777 


6 . 0 


2.4998 


4 . 7740 


1.5568 


3.2172 


7.0 


2 . 5000 


8.5235 


1.5568 


6.9967 


8.0 


2.5000 


11.9169 


1.5568 


10 . 3601 


9.0 


2.5000 


14 . 3046 


1.5568 


12.7478 


10. 0 


2.5000 


15.6965 


1.5568 


14.1397 


20.0 


2.5000 


16.9043 


1.5568 


15.3475 
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+ (K 28 + K 29 t +K 30 t 2 +K 31 t 3 )e Bt 



(23) 



The inverse Laplace transform of additional terms of 
Equation (22) results in a steady-state solution given by 
Equation (20). The region in which the solution can be 
applied is then the same as discussed in the previous 



h = 1.0, the solution represented by Equation (23) and the 
solution obtained considering the first three terms of 
Equation (23) are compared to a machine solution in Table 
IV. The small errors listed indicate that a highly accurate 
approximation was obtained. Table V illustrates the effect 
of increasing the coefficient of the non-linear term to ten. 
A less accurate approximation was obtained. The effect of 
increasing the magnitude of the forcing function was not 
illustrated; however for the parameter values used in 
Table IV if c had been increased to ten, a steady-state 
value of -1.4063 would have been obtained indicating that 
the solution limits would have been exceeded. The linear 
solution is seen to be a less accurate approximation. 

3 . Radical Imaginary 



2 

section. Using the values c = 1.0, a = 6.0, Wq = 4.0 and 



2 2 

The relationship between a and 4 wq is: 



„ 2 . 2 
4oJq > a 



(24) 
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TABLE IV 





a =6 . 0 } 


Wq= 4 . 0 , h=l . 0 , 


O 
1 — 1 

II 

O 




Time 


Linear 


Reversion 


Machine 


Error 


0.0 


0.0000 


0.0000 


0.0000 


0.0000 


0.5 


0.0533 


0.0533 


0.0533 


0.0000 


1.0 


0.1139 


0.1139 


0.1138 


0.0001 


1.5 


0.1570 


0.1569 


0.1568 


0.0001 


2.0 


0.1865 


0.1864 


0.1860 


0.0004 


2.5 


0.2066 


0.2064 


0.2058 


0.0006 


3.0 


0.2204 


0.2200 


0.2191 


0.0009 


3.5 


0.2298 


0.2292 


0.2281 


0.0011 


4.0 


0.2362 


0.2353 


0.2341 


0.0012 


5.0 


0.2436 


0.2420 


0.2408 


0.0012 


6.0 


0.2470 


0.2448 


0.2438 


0.0010 


7.0 


0.2486 


0.2459 


0.2452 


0.0007 


o 

oo 


0.2494 


0.2462 


0.2459 


0.0003 


9.0 


0.2497 


0.2463 


0.2460 


0.0003 


10.0 


0.2499 


0.2463 


0.2462 


0.0001 


11.0 


0.2499 


0.2462 


0.2462 


0.0000 


12.0 


0.2500 


0.2462 


0.2463 


- 0.0001 


13.0 


0.2500 


0.2461 


0.2463 


- 0.0002 
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TABLE V 





P 

It 

C T\ 
O 

V 0 


Wq = 4 . 0 , h= 10 . 0 , 


O 

l — 1 

II 

O 




Time 


Linear 


Reversion 


Machine 


Error 


0.0 


0.0000 


0.0000 


0.0000 


0.0000 


0.5 


0.0533 


0.0533 


0.0533 


0.0000 


1.0 


0.1139 


0.1139 


0.1136 


0.0003 


1.5 


0.1570 


0.1568 


0.1553 


0.0015 


2.0 


0.1865 


0.1857 


0.1821 


0.0036 


2.5 


0.2066 


0.2047 


0.1986 


0.0061 


3.0 


0.2204 


0.2165 


0.2086 


0.0079 


3-5 


0.2298 


0.2234 


0.2145 


0.0089 


4.0 


0.2362 


0.2269 


0.2179 


0.0090 


4.5 


0.2406 


0.2281 


0.2196 


0.0085 


5.0 


0.2436 


0.2278 


0.2210 


0.0068 


5.5 


0.2465 


0.2266 


0.2216 


0.0050 


6.0 


0.2470 


0.2250 


0.2220 


0.0030 


6.5 


0.2480 


0.2231 


0.2222 


0.0009 


7.0 


0.2486 


0.2213 


0.2223 


- 0.0010 


8.0 


0.2494 


0.2181 


0.2224 


- 0.0033 


9.0 


0.2497 


0.2157 


0.2225 


- 0.0068 


10.0 


0.2499 


0.2139 


0.2225 


- 0.0086 


12.0 


0.2500 


0.2120 


0.2225 


- 0.0109 


15.0 


0.2500 


0.2110 


0.2225 


- 0.0115 
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Let : 



1 

2 



(a 2 -4w 2 ) = jD. 



Then : 



A = | + jD (25) 

and 



B = 




( 26 ) 



Equation (15) is then expressed as: 



Y(s) = 



he 

Tj 



,, 2 3 
3h c^ 



s(s+A)(s+B) s(s+A) (s+B) s ( s+A) ^ ( s+B ) ^ 



(27) 



The inverse Laplace transform of the first two terms of 
Equation (27) yields: 



y(t) =K^ 0 +K 4l e At + K^ 0 e Bt + + (K^ + K^t + K^t 2 +K il7 t 3 ) e Afc 



L 42 



il 3 T ^44 T ^45' 



L 46' 



k 47 



+ (K 4 g +K 49 t +K 5 Q t 2 +K 51 t 3 )e Bt 



( 28 ) 



Since A and B are complex, the coefficients of the time- 
dependent portion of Equation (28) are complex. Equation 
(28) can then be rewritten as: 



y(t) = K 4o + (M 1 +jM 2 )e At +(M 1 -jM 2 )e Bt +K 43 

P P —At - 

+[(M 3 +jM 4 ) + (M 5 +jM g )t + (M ? +jMg)t^ + (M 9 +jM 1Q )t^]e 
+[(M 3 -jM 4 ) + (M 5 -jMg)t + (M 7 -jMg)t 2 +(M 9 -jM 10 )t 3 ]e _Bt 
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( 29 ) 



Combining terms Equation (29) becomes: 

_«t --t 

y (t ) = K^ 0 + 2M x e 2 cos(Dt) + 2M 2 e 2 sin(Dt) 

_« t 

+ (2M 0 +2M [ -t + 2M 7 t 2 + 2M n t 3 )e 2 cos(Dt) 

o b ! y 

+ (2M 1| + 2M 6 t + 2Mgt 2 + 2M 10 t 3 )e 2 sin(Dt) . (30) 

In this case as in the previous two cases considered, the 

steady-state solution is given by Equation (20). In Table 

VI, the reversion solution is compared to a machine solution 

for a = 1.0, Wq = 4.0, h = 1.0 and c = 1.0. In Table VII, 

the solutions are compared for the non-linear coefficient 

equal to ten. A comparison of the results listed in Tables 

VI and VII with the results obtained in the previous cases 

2 

for the same values of h and c, indicates that the 

reversion method yields a less accurate transient solution 
when the damping coefficient is small. In this case the 
linear solution is seen to be a better approximation over a 
portion of the transient solution. 

B. SINUSOIDAL FORCING FUNCTION 

Substitution of a constant-amplitude sinusoidal forcing 
function, c sin(iot), into the first two terms of Equation 
(15) and utilization of the trigonometric identity [Ref. 63 

sin 3 (wt) = ■^-sin(wt) - jj-sin(3wt) (31) 

yields : 
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TABLE VI 



a=1.0, Wq= 4.0, h=1.0, c=1.0 



Time 


Linear 


Reversion 


Machine 


Error 


0.0 


0.0000 


0.0000 


0.0000 


0.0000 


0.5 


0.0982 


0.0982 


0.0982 


0.0000 


1.0 


0.2677 


0.2677 


0.2672 


0.0005 


1.5 


0.3576 


0.3574 


0.3534 


0.0040 


2.0 


0.3343 


0.3330 


0.3235 


0.0095 


2.5 


0.2591 


0.2550 


0.2463 


0.0087 


3.0 


0.2069 


0.1985 


0.1994 


-0.0009 


3.5 


0.2064 


0.1951 


0.2060 


-0.0109 


4.0 


0.2377 


0.2279 


0.2402 


-0.0123 


4.5 


0.2655 


0.2620 


0.2655 


-0.0035 


5.0 


0.2712 


0.2743 


0.2663 


0.0080 


5.5 


0.2593 


0.2642 


0.2516 


0.0126 


6.0 


0.2453 


0.2453 


0.2387 


0.0066 


6.5 


0.2403 


0.2324 


0.2367 


-0.0043 


7.0 


0.2442 


0.2320 


0.2426 


-0.0106 


8.0 


0.2542 


0.2510 


0.2508 


0.0002 


9.0 


0.2503 


0.2523 


0.2454 


0.0069 


10.0 


0.2483 


0.2409 


0.2450 


-0.0041 


50.0 


0.2500 


0.2461 


0.2463 


-0.0002 
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TABLE VII 





a=l . 0 , 


O 

O 
1— 1 

II 

O 

II 

C\J 0 
3 


O 
1 — 1 

II 

a 




Time 


Linear 


Reversion 


Machine 


Error 


0.0 


0.0000 


0.0000 


0.0000 


0.0000 


0.5 


0.0982 


0.0982 


0.0982 


0.0000 


1.0 


0.2677 


0.2675 


0.2626 


0.0049 


1.5 


0.3576 


0.3554 


0.3188 


0.0366 


2.0 


0.3343 


0.3211 


0.2507 


0.0704 


2.5 


0.2591 


0.2173 


0.1799 


0.0374 


3.0 


0.2069 


0.1230 


0.1761 


- 0.0531 


3.5 


0.2064 


0.0934 


0.2164 


- 0.1230 


4.0 


0.2377 


0.1400 


0.2463 


- 0.1063 


4.5 


0.2655 


0.2298 


0.2410 


- 0.0112 


5.0 


0.2712 


0.3026 


0.2195 


0.0831 


5.5 


0.2593 


0.3085 


0.2093 


0.0992 


6.0 


0.2453 


0.2451 


0.2156 


0.0295 


6.5 


0.2403 


0.1619 


0.2260 


- 0.0641 


7.0 


0.2442 


0.1222 


0.2289 


- 0.1067 


8.0 


0.2542 


0.2225 


0.2198 


0.0027 


9.0 


0.2503 


0.2703 


0.2222 


0.0481 


10.0 


0.2483 


0.1745 


0.2234 


- 0.0489 


50.0 


0.2500 


0.2110 


0.2225 


- 0.0115 
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